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Abstract A theoretical analysis of Pancharatnam and Berry phases is 
made for biphoton three-level systems, which are produced via frequency 
degenerate co-linear spontaneous parametric down conversion (SPDC). The 
general theory of Pancharatnam phases is discussed with a special emphasis 
on geodesic 'curves'in Hilbert space. Explicit expressions for Pancharatnam, 
dynamical and geometrical phases are derived for the transformations pro- 
duced by linear phase-converters. The problem of gauge invariance is treated 
along all the article. 
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1. Introduction 



An interesting development in quantum mechanics has been made by 
Berry [1] who has shown that when a quantum-mechanical state is devel- 
oping around a closed circuit it obtains a geometrical phase in addition to 
the well known dynamical phase. Aharonov and Anandan [2] have general- 
ized the use of geometric phase to non-adiabatic time developments. Certain 
ideas from optics by Pancharatnam [3] have been used for defining also phase 
changes for partial and non unitary cycles [4,5]. The discovery of a topolog- 
ical phase in quantum mechanics has led to a unified approach to various 
topological phenomena in physics, both at the quantum and classical levels. 
The connections between the geometric phase and fibre-bundle theories [6-10] 
have been analyzed. Following Berry's paper many experiments have been 
developed for measuring the geometric part of the phase change, in a variety 
of contexts, which are related to Lie group theories. Most observations of 
Berry's phases are related to SU(2) and Lorentz groups, and there is an enor- 
mous amount of literature on such systems. In a previous paper [11] Berry's 
phases and interference effects for three- level atoms, related to the SU(3) al- 
gebra, have been analyzed, and the comparison between this work and those 
of other authors have been discussed. In the present article we would like to 
develop the theory of Pancharatnam and Berry phases for three-level optical 
systems related to the SU(3) algebra. More specifically the theory of Pan- 
charatnam phase is developed here for a three-level photonic system which is 
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realized by the polarization states of the single mode biphoton field [12-15]. 
A general review for the unitary evolution of polarized light governed by the 
SU(2) group and its relation to topological effects has been published [16]. 
We extend such treatment to three-level photonic system. 

There are close connections between the two concepts of Pancharatnam 
and Berry phases. Since we are treating an optical system we find it more 
useful to base our analysis on the concept of Pancharatnam phase, but we 
would also like to show the connection of our analysis with the concept of 
Berry phase. For the basic understanding of our analysis for a specific system 
we find it useful to review shortly the basic theories which are used in our 
derivations. 

The paper is arranged as follows: In Section 2 we show the relation be- 
tween Berry phase and topology. In Section 3 we describe various properties 
of the Pancharatnam phase and its relation to the Berry phase. The main 
part of the paper is given in Section 4, where we analyze possible interfer- 
ence effects which can be obtained in single mode biphoton fields [12-15]. In 
Section 5 we summarize our results. 

2. Berry Phase Related to Topology 

When a quantum-mechanical (QM) state is evolving adiabatically in time 
and transported around a circuit as an eigenstate of a Hamiltonian with 
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slowly varying parameters, it acquires a geometrical phase factor in addition 
to the dynamical phase [1]. The geometrical phase factor is given by 



7n (T) = i j dt(n, R(t) | d/dt | n, R(t)) = i j> dR(t) ■ (n, R(t) | | n, R(t)) (1) 

for a Hermitian Hamiltonian H(R(t)) with parameters Ri(t) which are slowly 
varying along a closed curve C, in parameter space in time T. In the adi- 
abatic approximation, it is assumed that | n, R(t)) is an eigenstate of the 
Hamiltonian: 

H(R(t)) | n, R(t)) = E n R(t) | n, R(t)) (2) 

The phase factor acquired by the state | ^ n (R,T)) consists of two parts: 

| V n (R, T)> = exp(z 7n (T) exp(-^ F E n (r)dr) \ n, R(t)) (3) 

Jo 

where the first and second exponents on the right side of Eq. (3) represent, 
respectively, the adiabatic Berry phase and the dynamical phase. The term 

A M {n,t)=i(n,R(t) \V n \n,R(t)) (4) 

is referred to as a vector potential. As in electromagnetic theory the vec- 
tor potential is defined up to a gauge transformation. By performing the 
transformation 

| n , R(t))' = exp{ia n {R))\n, R(t)) (5) 
we induce a "gauge transformation" on Ag(n, t): 

A' n (n,t) = A A (n,t)-V A (a n (R)) (6) 
4 



/,From the definition of Berry phase given by Eq. (1) it is clear that the 
Berry phase is gauge invariant for a closed circuit. 

Considering the state vector | ^f n (R,t)), Aharonov and Anandan [2] re- 
moved its dynamical phase factor, by defining a new state: 

|0„CR,f) = ex P[* t h n (R,t')dt']\y(n,R(t)) (7) 
Jo 

where 

— * 

\4> n (R,t')) satisfies the equation: 

id/dt\<j> n (R,t)) = [H(t) - h n (R,t)}\cj> n (R,t)) (9) 

Multiplication of this equation on the left by ((f) n (R,t)\ yields the 'parallel 
transport' law: 

Im((f> n (R,t)\d/dt\(f> n (R,t)) = (10) 

Following fibre-bundle theories [6-10] we refer to Eq. (10) as the 'connec- 
tion'. The 'total space' N, in our analysis, includes all normalized Schrodinger 
wavef unctions. As is well known it is possible to multiply the wavefunction 
by a phase factor exp(i#), without changing the physical properties of the 
quantum state. All wavefunctions which differ only by a phase factor are 
considered in QM as one ray. The 'base space' R, in the present analy- 
sis, includes all rays of Schrodinger wavefunctions. There is a projection 
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7r : N — > i? from the total space to its base, i.e. from each wavefunction 
onto its ray. The 'fibre space' F is related to our analysis by including all 
wavefunctions which differ only by a phase factor as one 'fibre'. 

Let us assume, for simplicity, that an element of the 'bundle' \(f> n (R,t)) 
can be written 

\cf> n (R,t)) = n\z n (R,t)) (11) 

where \z n (R,t)) is representative of a certain ray and Q is the phase factor, 
Q = exp(i0). On the right side of Eq. (11) Q represents the 'vertical' 

— * 

part, which is along the 'fibre', while \Z n (R,t)) represents the 'horizontal' 
part along the 'base'. Substituting Eq. (11) into the connection, Eq. (10), 
performing the derivative and integrating this equation over a closed circuit 
during time T we have two contributions which are equal in magnitude and 
opposite in sign, so that we obtain: 

Arg(fi(T)) - arg(ft(0)) = ln{T) =ijdR- {(f> n (R, t)\V a \(j> n (R(t)) (12) 

The geometrical phase factor is then given in Eq. (12) which is the same 
relation as that of Eq. (1), but without any adiabatic approximation. One 
finds according to the above analysis that although the quantum circuit is 
closed in the 'base space' (\z n (T)) = |z n (0))) it is open in the 'total space', 
as \<f) n (T)) is different from |</> n (0)), and this difference is given by the Berry 
phase. 
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3. Pancharatnam Phase and its Relation to Topological Effects 
a. Definition of Pancharatnam Phase 

Consider two normalized nonorthogonal Hilbert states \A) and \B), and 
assume further that \A) is exposed to £7(1) shift [17]. The resulting 
interference pattern is determined by 

/ = \e^\A) + \B)\ 2 = 2 + 2\(A\B)\ cos[0 - aig{A\B)] (13) 

where its maximum is obtained at the Pancharatnam relative phase 0o = 
arg{A\B). This phase is reduced to the £7(1) case whenever \B) = e %e \A) 
as it yields wg(A\B) = 9. In the original treatments of the Berry phase 
[l,2]one considers a quantal system evolving around a closed circuit, from an 
initial wavefunction \A) to a final wavefunction \B) where \B) is obtained 
from \A) by a cyclic evolution, i.e., by multiplication with a U(l) phase fac- 
tor. Although the initial phase of the quantal state (its 'fibre') is defined 
arbitrarily, the phase difference between the state \A) and \B) is well defined 
and can be observed by interferometric methods. Pancharatnam provided 
the physical insight [17] that in a noncyclic evolution if the final wave \B) 
is superimposed on the initial state \A) only the component (A\B)\A) along 
| A) interferes with \A). All other components of \B) which are orthogonal 
to | A) merely add to the intensity, since their cross terms with \A) vanish. 
The Pancharatnam 'connection' defines the phase between \B) and \A) as 
4> Q = axg{A\B). The interference amplitude differs from unity if the 

evolution is non-cyclic [17,18], and this difference leads to reduction of the 
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visibility in the interference pattern. The Pancharatnam phase is indeter- 
minate if \B) is orthogonal to \A). Although Pancharatnam phase can be 
defined also for mixed states [19-21] we restrict the analysis of the present 
paper to pure states. The use of Berry and Pancharatnam phases to non- 
cyclic evolution raises the problem of gauge invariance which is treated in 
the following paragraph. 

b. Pancharatnam Geometric and Dynamical Phases 

An extensive discussion of the mathematical and physical properties of 
the Pancharatnam phase has been presented in the articles of Mukanda and 
Simon [22]. We describe here only a few fundamental properties which are 
needed for understanding the analysis given in Section 4 for a specific system. 

Consider a one dimensional or one parameter 'curve' consisting of a family 
of wave vectors |^(s)) which is changing continuously as a function of s. s 
is any parameter by which \l/(s) is continuously changing (including time as 
a special case). Assuming that |^/(s)) is a unit vector for all |^(s)) we get 

Re(V(s)\V(s)) = 0; (*(s)|*(s)) = i Im(V(s)\if(s)) (14) 

where the derivatives are according to the parameter s. One can change the 
wavevectors |\&(s)) by a "gauge transformation": 

|tt(s))^|tt'(s))=e ia «|tt(s)); se[ Sl ,s 2 ] (15) 



8 



where a(s) is a real smooth function of s in the interval se[si, S2]. For cyclic 
transformation described in the previous section the Berry phase is indepen- 
dent of the gauge transformation. However, for noncyclic transformation we 
get 

Im(y'(s)\&(s)} = Jm(tf (s), *(s)) + d(s) (16) 

^From the above equations we find that we can construct a functional in 
Hilbert space which is gauge invariant: 

arg(*'(si),*'(s 2 )) -Im H ds(*'(s)|*'(s)> = (17) 

Jsi 

arg(\E'(s 1 ), ^(s 2 )) — Im / ds (^(s)|^(s)) = gauge invariant 

This property can be interpreted according to the topological description 
in the previous section. Here the wavefunction ^(s)) includes its 'fibre' 
i.e. |^(s)) belongs to the total space N where it includes all wavefunctions 
which differ in a phase factor exp(i#) [changing this phase factor by the gauge 
transformation is considered in topology as moving along the 'fibre']. On the 
other hand, the functional appearing in Eq. (17), which is gauge invariant 
represents a projection from the total space to its basis and is defined as 
the "geometric phase". We find according to Eq. (17) that by substracting 
the "dynamical phase", which is defined as Im ds{if(s)\4f(s)), from the 
Pancharatnam phase, which is defined here as arg( , I'(si), \l/(s2)), we get the 
geometric phase which is gauge invariant. The interesting point here is that 
while Berry phase analysis, as given in Section 2, is well defined (gauge in- 
variant) only for cyclic evolution, the geometric phase is well defined also for 
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noncyclic evolution (also for mixed states [19-21]). Eq. 17 can be reformu- 
lated as [22]: 



where <p p (N) and (pd y n(N) are defined over the total space N and each of 
them is gauge dependent while <fi g (C) is the geometric phase which is defined 
over the basis space of rays and is gauge invariant. Eqs. (17) and (18) have 
been restricted to a dependence on only one parameter s, but such restriction 
is valid for the optical system which will be treated in Section 4. For a change 
of the wavefunction in an interval [s±, s 2 ] one finds that the Pancharatnam 
phase is fixed only by the initial and final wavefunctions |\&(si)) and |\P(s 2 )), 
respectively. In order to find the geometric phase one should find a way to 
subtract the dynamical phase which is defined at each point s and is obtained 
by its integration over the interval [s±, s 2 ]. 

One can measure the Pancharatnam phase by applying Eq. (13), but 
it might be quite difficult to subtract experimentally from it the dynamical 
phase, although it can be done by using a theoretical calculation. In consid- 
ering consecutive wavefunctions 'curves' one should take into account that 
the geometric and Pancharatram phases do not have the additive property 
[19]. In comparison, the dynamical phase of the total path can be obtained 
by adding the dynamical phases of the 'curves' from which the total path 
is composed. Especially interesting are the geodesic 'curves' that extremize 



4> g {C) = <j> p (N) - dyn (AO; M N ) = arg{(*( Sl )|M>(s 2 )}}; (18) 
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the distance with respect to a certain metric. In the next paragraph we in- 
troduce a metric for the present Hilbert space and show how such extremum 
can be obtained. The geometric phase which is obtained by subtracting the 
dynamical phase from the Pancharatnam phase can have global properties 
while the geodesic 'curve' has also local differential properties. These prop- 
erties are obtained by extremizing a certain functional, leading to the path 
with the shortest distance between two points in Hilbert space. 

c. Hilbert Space Metric and Geodesies 

The distance between two quantum states |^(s)) and ^(s')), representing 
a quantum system which is developed according to parameter s, can be given 
as [23,24] 



As mentioned previously s might represent not only the time t but also 
any other parameter s by which the wavefunction is continuously developed. 
The use of Eq. (19) is reasonable since if the ray |^(s)) is orthogonal to the 
ray |^(s')) then the distance between the two states is equal to 1, while if 
|^(s')) 1S equal to |^(s)) the distance vanishes. Additional explanation to 
this definition will be given also later, but let us see first the mathematical 
derivations which are obtained from this equation. 
Assuming two close wavefunctions |^(s)) and + ds)), then we get: 




(19) 
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By using the second derivative of (^(s)|^(s)( we get 




(21) 



By using Eqs. (19-21) one gets, after some algebra [23,24]: 

(^f) ={l-(^(s)\^(s + ds))(^(s + ds)\^(s))}/ds 2 (22) 




The derivation of Eq. (22), which has been made by using the definition 
(19), can be obtained also in a different way. The horizontal component of the 
tangent vector d/ds\^(s)) is given by d/ds\^(s))-{^(s)\^(s))\^(s)). Here 
we have subtracted from the derivative of the wavefunction its movement 
along the 'fibre' [25], since it does not change the basis of the wavefunction 
i.e., its ray. The norm of the above vector is given after a straightforward 
algebra as: 



which is equivalent to Eq. (22) up to a multiplication by factor 2, which can 
be inserted arbitrarily in the definition (19). This equation gives a physical 
insight for using Eq. (19). 

In deriving Eq. (22) one neglects terms which are of order ds 3 or higher. 
Given the continuous 'curve' C of the wavefunction ^f(s) varying continuously 
as a function of the parameter s from si to s 2 one gets the functional [22,26]: 





(23) 
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The geodesic curve is obtained by extremizing this functional. Mukanda and 
Simon [22] have shown that any 'curve' composed of a normalized wavefunc- 
tion which is changing continuously as a function of s and which obeys the 
equations: 

|*( a )> = -(*( a )|*( a )>|*( a )>, (24) 



(*(a)|*(a)> = 0, (25) 
is a geodesic curve. The path [22, 26-28] 

|*(a)> = |A> oos(a) + l l ^~ lA){AlB) \sms (26) 

I y/l~\(B\A)\ 2 J 

where \A) and \B) are the initial and final wavefunctions is a geodesic curve 
as it fulfills Eqs. (24) and (25). 



The use of the geodesic equations (26) for noncyclic £77(2) evolution has 
been discussed [27]. In a recent interesting article [28] this equation has been 
applied for discussing possible geometric phase measurement of three-level 
systems in interferometry. According to Eq. (25) \^f) is 'parallel transported' 
[see Eq. (10)] as it leads to vanishing of the dynamical phase at each point s. 
The 'horizontal' property given by Eq. (25) can be destroyed, but not the 
geodesic property, by multiplying the function |^(s)) by e ia ^ where a(s) 
is any continuous function of s. Such gauge transformation does not destroy 
the geodesic property as a geodesic is gauge invariant (any property of the 
basis of the wavefunctions i.e., of the rays is gauge invariant). 
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d. The Vertex Theorem 



Let us assume that we have N — 1 consecutive wavefunctions 'curves' 
where the r'th 'curve' is described by the transition \l/ r — > \l/ r+1 and the total 
path is obtained by summation over r from 1 to N — 1. Using Eq. (18) we 
get: 

N-l 

4> a {C) = arg(* 1 |* iV ) - £ dyn (M/ r - * r+1 ) (27) 

r=l 

Although the Pancharatram and the dynamical phase can be gauge depen- 
dent their difference given as <f> g (C) is gauge independent. 

For cases in which the geometric phase vanish in all the transitions \EV — > 
^ r +i we get: 

N-l 

<P g (C) = arg(^ 1 |^ 7V ) - Yl arg(^ r |^ r+ i) = 

r=l 

-axg{(*jv|*i)(*i|* 2 )(*2|*3)-(*jv-i|*jv>} (28) 
Eq. (28) can be considered as the general "vertex theorem" . 

Mukanda and Simon [22] have proved, by developing certain mathemati- 
cal procedures, that any two normalized vectors in ray space can be connected 
by a geodesic curve. Therefore one can apply Eq. (28) for the special case 
in which each of the transitions \l/ r — > ^ r +i is a geodesic. But in deriving 
Eq. (28) we have assumed only that the geometric phase in the "global" 
transition \l/ r — > ^ r+ i vanish and this does not necessarily imply that this 
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transition is geodesic. However, if we assume that the geometric phase vanish 
continuously for each of the differential transition — > \I/(s + e) (e — > 0) 
in the interval se[sis 2 ] then the 'curve' is geodesic and Eq. (28) gets the form 
[20]: 

UC) = 

- axg{(*(s 2 )|*(s 1 )>(*(s 1 )|*(s 1 + e))(*(sx + e)|*( Sl + 2e))... 
(*( Sl + (iV-l)e)|*( S2 ))} (29) 

where (N — l)e = s 2 — s 1 , (N — 1) — > oo , e — > 0. 

Due to the fact that 4> g (C) is gauge invariant one can use "parallel lift" 
by which (\[>(s)|\[>(s)) = and then the vanishing of the geometric phase in 
each of the intervals (^(s), *f?(s + e)) implies (s) ^(s + e)) = 1. Then we 
get the simple expression 

9 (C) = -arg(M/( S2 )|M/( Sl )> . (30) 

For such cases one may measure 4> g (C) as a relative phase shift in the inter- 
ference pattern by applying Eq. (13), assuming \A) = \ty(si)), \B) = \^(s2)). 
However, if the development is not given by 'parallel transport' one has to 
subtract the dynamical phase on the right side of Eq. (30). 



15 



4. Pancharatnam Phase for Polarized Biphotons 
a. Descriptions of the Biphoton States 

We consider a quantum system formed by two correlated photons - a 
biphoton, emitted via frequency degenerate colinear spontaneous parametric 
down conversion (SPDC) [12]. We assume that the biphoton state can be 
described as 

|*) = c 1 |2,0) + c 2 |l,l) + c 3 |0,2) (31) 

where (N x , N y ) means a state with N x photons in the horizontal (x) polar- 
ization mode and N y photons in the vertical (y) polarization mode, with 
N x + N y = 2. The states |2, 0) and |0, 2) are generated via type-I SPDC and 
the state |1, 1) via type-II SPDC. Arbitrary transformations of the polariza- 
tion vectors (ci,c 2 ,c 3 ) [| c x | 2 + c 2 | 2 + |c 3 | 2 = 1] are given by unitary 3x3 
matrix G, where G^G = I, det G — 1, which form a three-dimensional rep- 
resentation of the SU(3) group. 

By passing from the basis |2, 0), |1, 1), |0, 2) to the basis 

= |,. > = M^M ;|So) = |u) (32) 

one obtains three states that can be transformed into one another by means 
of only phase-plates. It has been suggested to use polarized biphotons as 
ternary analogs of two-state quantum systems (qubits) [12,15]. Our aim in 
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the present paper is, however, different. We would like to study here pos- 
sible interference effects, by the use of biphotons which can be related to 
Pancharatnam and Berry phases in a three-level photonic system. The effect 
of the loss-free polarization converters (phase-plates) on a biphoton state 
given by Eq. (31) has been described by Burlakov and Klyshko [14]. For 
our purpose the effect of the phase-plates is described by a 3 x 3 unitary 
matrix which operates directly on the basis of states given by Eq. (32). This 
transformation matrix and its properties are described in the next paragraph. 



b. Transformations of Biphotons by Phase-Plates 



The effect of polarization converters on a biphoton given by Eq. (31) as 
(ci,c 2 ,c 3 ) is described by the transformation matrix [14]: 



G = 



( 



V 



t 2 

-V2tr* 



V2tr 
\t\ 2 -\r\ 2 
-V2t*r* 



r 2 \ 
V2t*r 

i*2 



t* 



(33) 



Here t and r are amplitude transmission and reflection coefficients of a given 
converter. We assume a linear phase-plate with optical thickness S and orien- 
tation x relative to the horizontal direction x which corresponds to the trans- 
formation t = cos 5 + isin<5cos(2x), r = i sin 5sin(2x). For a quarter- wave 
plate we have 5 = 7r/4 and then t — (1 + ?cos(2x))/v^2, r = i sin(2x))/v / 2. 
A half-wave plate gives t — i cos(2x), r — i sin(2%). 
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The basis of states |2, 0), |11) and |0,2) can be transformed to the basis 
of states |^+), |^_) and |^ ) as 



(34) 







( 2,0 \ 


^_ 


= A 


1,1 


V ^0 J 




V 0,2 J 



where the matrix A is given by 



A 



It* 75 \ 

J_ n 1 - 

V2 U V2 



(35) 



The normalized biphoton state can be expressed in the new basis by 



(36) 



where the effect of the polarization converters on the biphoton of Eq. (36) 
given as (d±, d 2 , d 3 ) is described by the transformation matrix: 



Q = AG A' 1 = 



I t 2 +r 2 +r* 2 +t* 
2 



2 j_, r 2 _ r .*2 _ f *2 t 2_ r 2_ r *2_|_ t *2 



t z +r z -r* 2 -t 



tr — t*r* \ 
tr + t*r* 



\ -tr*+t*r -tr*-t*r \t\ 2 - \r\ 2 ) 



(37) 



Assuming a linear phase-plate with optical thickness S and orientation \ the 
matrix Q gets a more explicit form as: 

/ cos(2<5) % sin(2<5) cos(2x) i sin 5 sin(2x) \ 



Q 



i sin(25) cos(2x) cos 2 5 — sin 2 5 cos(4%) — sin(4x) sin 2 S 
V i sin(25) sin(2x) — sin(4x) sin 2 (5) cos 2 5 + sin 2 5 cos(4x) / 



(38) 
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An interesting property of the matrix of Eqs. (37) and (38) is that its diago- 
nal matrix elements are real. The relations between the matrix Q operating 
on the biphoton state (di,d 2 ,d 3 ) and Pancharatnam and Berry phases will 
be the subjects of the following discussions. 

c. Biphoton Dynamical and Pancharatnam Phases 

Following the general analysis presented in Section 3 the Pancharatnam 
phase is developed as a function of the parameter s and this is analogous to 
the development in time of Berry phase by Schrodinger equation. In order to 
apply the general theory to phase converters we identify the parameter s as 
the parameter <5(or later as 25) in the transformation matrix Q given by Eq. 
(38). The input state vector is defined as (d 1 ,d 2 ,d 3 ) while the output state 
is obtained by multiplying this input state by the matrix Q. We assume 
that the initial state is developed continuously by changing in the matrix 
Q the parameter 5 from its initial value zero (giving a unit matrix), to its 
final value 5q, which corresponds to the optical depth of the converter. By 
straightforward calculations we find that the linear phase converter leads to 
a dynamical phase shift given by: 



/ di \ 




(39) 



V 



d 3 



J 




d5{2icos(2x)[did* 2 + d\d 2 ] + 2i sm.{2x)[d\d^ + d\d 3 }} 
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The result obtained in Eq. (39) is quite simple showing that the integrand 
in this integral is constant, depending only on the initial state (d±, di, ds) 
and the orientation x which has a fixed value for a certain converter. The 
dynamical phase is obtained by multiplying this integrand by 5 . 



'Parallel transport' of the biphoton state vector is obtained by the re- 
quirement that the integrand of Eq. (39) vanishes. 'Parallel transport' of 
the biphoton gives, however, trivial results since it is easy to show that un- 
der this condition both the Pancharatnam and geometric phases vanish. This 
conclusion can be obtained by calculating the imaginary value of (^ml^out) 
where \^m) is the initial biphoton state vector while |^ ut) is obtained by 
multiplying this initial state by the transformation matrix Q. We get: 

/ d ± \ 

Im(^ m \q out ) = Im{(d* 1: d* 2 , d\)Q d 2 } = (40) 

\ds J 

sin2<5{cos2x(d*d 2 + d* 2 di) + sm{2x){d* 1 d z + did\)} 

Therefore, vanishing of the integrand in Eq. (39), which is the condition for 
'parallel transport' implies vanishing also of the Pancharatnam phase. Cases 
for which the total path is composed by a continuum of differential 'curves' 
elements, where for each of them the geometric phase vanishes while the geo- 
metric phase for the total path might be different from zero can be defined as 
geodesic curves in its most general meaning. For an infinitesimal evolution of 
the biphoton state from an initial state (d±, d 2 , d 3 ), obtained by using addi- 
tional thickness d5, the change in geometric phase vanishes. This conclusion 
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can be obtained by using in Eq. (40) the approximation sin (25) ~ 25 and 
compare the Pancharatnam phase change with that obtained for the dynam- 
ical phase change given by Eq. (39). Therefore, the general evolution of the 
Biphoton state by phase converters can be considered as geodesic curves in a 
general sense. More basic and restricted definitions of geodesic 'curves' will 
be given later in the article. 

There are other interesting effects which can be observed by linear phase 
converters: 

1) By calculating the eigenvalues and eigenvectors of the transformation 
matrix Q one can describe cyclic transformations by which the initial state 
vector is multiplied by U(l) phase factor. 

2) Special geodesic developments can be obtained for biphoton states under 
special conditions, which can lead to interesting interference effects. 

3) Under general conditions, the geometric phase can be calculated by using 
the general expressions of Eq. (18). 

d. Eigenvalues and Eigenvectors of the Transformation Matrix 

Since the transformation matrix Q is unitary its eigenvalues are given 
by complex numbers with a unit absolute value. By calculating the eigen- 
vectors of the matrix Q one finds the initial biphoton states for which the 
transformation Q produces cyclic transformations i.e., the initial state vec- 
tors are multiplied by these complex numbers. The general calculation of 
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the eigenvalues of the matrix Q for arbitrary optical depth 5 leads to cubic 
equations. For simplicity, we demonstrate the general procedure by applying 
it to relatively simple cases. 

For a quarter- wave plate with 26 = 7r/2 the determinant for the eigenval- 
ues A is given by 

-A 3 + A 2 - A + 1 = (41) 

The roots of Eq. (41) are given by A = ±i, 1 (independent of the orientation 
x)- For the eigenvalue A = % the normalized eigenvector is given by 

1 cos(2x) sin(2x) , 

di = ^, 4. = -^. d * = -jT (42) 

Here we have assumed that d\ is real and this can be done since the phase 
of the initial biphoton state is arbitrary. However, the phase factor e llT l 2 by 
which the initial state is multiplied for getting the final state vector is well 
defined and can be observed in interference experiments. The dynamical 
phase can be calculated for the initial state, given by Eq. (42) with 25 = n/2, 
by Eq. (39) as: 

/" V4 d5{cos 2 (2 X ) + 2% sin 2 (2x)} = in/2 (43) 
J o 

So we find that all the Pancharatnam phase is a dynamical phase. Similar 
calculations can be given for the roots A = — i and 1. 



For a half-wave plate the straightforward calculations show that the roots 
are given by = ±1. This result is obvious since the half- wave plate can be 
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considered as the operation of two consecutive quarter-wave plates with the 
same orientation x-> so that the roots for the half-wave plate are obtained 
by the squares of the roots of the quarter-wave plates. One finds that the 
roots for phase-plates with arbitrary optical depth 5 will be given by e ±t2S , 1. 
One can also calculate the eigenvalues and the eigenvectors for multiplica- 
tion of two (or more) phase plates with different orientations x- Our interest, 
however, is to calculate geometric phases for noncyclic developments and for 
this purpose geodesic developments are described in the next paragraph, for 
special cases. 



e. Geodesic Development of Biphoton States 



Let us assume a phase converter with orientation given by cos (2%) = 1. 
We assume also that this phase converter transform a special initial biphoton 
state for which d 3 = (|<ii| 2 + \d 2 \ 2 = 1). We can easily prove that this 
biphoton state is developed as a function of 5 along a geodesic 'curve', as the 
geodesic equations (24) is obeyed. Defining s = 25 we obtain 



- sin s —i cos s \ / -sins % cos s 
-i cos s — sin s I \ i cos s — sin s 




(44) 



and as |^(s)) = — |^(s)), equation (24) is fulfilled. Although this bipho- 
ton does not obey the 'parallel transport' of Eq. (25) we can use a gauge 
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transformation of |^(s)) which will lead to a 'parallel transport'. This gauge 
transformation can be given according to Eq. (39) by 

|*(s)) -> |*'(s)) = |*(s)) exp{-is[di<£ + d*d 2 ]} (45) 

The conclusion from this transformation is that for the present biphoton state 
|^(s)) the geodesic property is preserved although the 'horizontal' property 
is destroyed, [see the discussions after Eq. (26)]. Since the 'horizontal' prop- 
erty is not preserved one can use Eq. (30) but the dynamical phase should 
be subtracted from the right side of this equation. One should therefore no- 
tice that the geodesic 'curve' described here is basically different from that 
suggested in Ref. [28]. 

For the present geodesic 'curve' we find that 9 = arg(\I/(si)|\I/(s2)) is given 
as a function of the optical depth 25 = s by 

tan 9 = tan(s) (d\d 2 + d*^) (46) 

and the geometric phase is given by 

<f> g {C) = 9- s{d\d 2 + d* 2 d ± ) (47) 

According to Eq. (46) 9 = s = 25 only for the special cases for which 
d\d 2 + d\d\ = 1, and only under this condition the geometric phase vanishes. 
For more general cases we find that when s is changing from 7r/2 — e (e — > 0) 
to 7r/2 + e there is a jump in the geometric phase of ir. Such phase jumps 
can be observed [16]. 
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In a similar way to the above analysis one can obtain another geodesic 
'curve' by using a phase converter with orientation given by sin(2x) = 1 and 
an initial state for which d 2 = (\di\ 2 + \d 3 \ 2 = 1). The present geodesic 
'curves' are obtained by the development of only two levels chosen from the 
three-level system. 

In the next paragraph we demonstrate calculations of the geometric phase 
for a real three-level system. We also explain the difference between the gen- 
eral biphotons transformation which we defined as geodesic 'curves' in the 
more general sense and the more restricted definition given by Eq. (24). 

f. Geometric Phase Obtained by Linear Converters 

By the derivations of Eqs. (39) and (40) we implied that the linear 
converters produce 'curves' which are geodesic in a general sense. Since the 
transformation of the biphoton state by the matrix Q of Eq. (38) does not 
obey the geodesic Eq. (24) we need to justify our assumption also from 
the mathematical point of view. The geodesic equation (24) is basically a 
harmonic oscillator equation for the wavefunction |\l/(s)). We find according 
to Eqs. (39) and (40) that only the imaginary elements of the matrix Q 
contribute to Pancharatnam and dynamical phases. The imaginary part of 
the Q matrix fulfills the equation 

Im{Q + 4Q) = (48) 
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which is basically an harmonic oscillator equation explaining the geodesic 
property of this transformation in its general sense. 



Let us demonstrate the use of Eqs. (18) for two consecutive transfor- 
mations operating on an initial state (d±, d 2 , d 3 ), first by using quarter-wave 
plate with an orientation x — and second by using half-wave plate with 
orientation x — 7r /4- We define the Pancharatnam phase in the first, second 
and total transformations as #p„ n , 9pi n , and 0pl n , respectively. Dynamical 
phases in the first, second and total transformations are defined as 9^ n , 9 d 2 y \ 

(3) 

and dyn , respectively. By straightforward calculations we get: 



3 (i) _ d\d 2 + d* 2 di _ a ( 2 ) _ a (tot) 

N 2 + \d 2 \ 2 + Wy/2 



tan 6 Pan — ■ - - ■ 1 i - - ■ 2 t - - j= ; 6 Pan — 6 Pan — ; 

eg, = 2n{{d\d z + d,di) + t(d 2 d; - did,)} ■ 

ft = K(d\d 2 + d* 2 di) ; 



dyn 

a(tot) _ fl (l) q{2) (aq , 
Vdyn - Vdyn + ® dyn I 49 J 

We find that the total Pancharatnam phase vanishes and therefore in the 
present case the total geometric phase is equal to minus the total dynamical 
phase. 



5. Summary and Discussion 



The present paper has analyzed the transformations of biphoton states 
which can be obtained by the use of linear phase converters in relation to Pan- 
charatnam and Berry phases. While most of the previous works have treated 
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two-level systems we analyze here a special three-level photonic system. In 
order to understand the present analysis the basic concept of Pancharatnam 
phase is explained in relation to topological effects. The fundamental prop- 
erties of Pancharatnam and Berry phases are reviewed for the purpose of 
using them in the analysis of a specific system. 

In Section 2 the relation between Berry phase and topology is explained. 
The phase obtained by an atomic system developing according to Schrodinger 
equation can be separated into a dynamical and geometric phase. This topo- 
logical separation is common for Berry and Pancharatnam phase effects in 
both atomic and photonic systems. The Berry phase is usually calculated 
for a closed circuit for which the calculation of phase is gauge invariant. An 
important development has been made by the use of Pancharatnam phase 
which gives, after the subtraction of the dynamical phase, a geometric phase 
which is gauge invariant. 

In Section 3 the main properties of the Pancharatnam phase are treated. 
The basic equation for obtaining the geometric phase in any atomic or pho- 
tonic system is given by Eq. (18). A special emphasis is made in the present 
work on geodesic 'curves' which are related to Hilbert space metric. The ba- 
sic equations for geodesies and for 'parallel transport' are given in Eqs. (24) 
and (25) respectively. A vertex theorem is developed for 'parallel transport' 
of a Hilbert state along a geodesic 'curve'. 
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The main results of the present work are given in Section 4 for analyzing 
Pancharatnam phases obtained by the transformations of biphotons using 
phase-plates. We have used the basis of states \& + , and \l/ given by Eq. 
(34) and the general transformation of these states by the phase-plates is 
given by Eq. (38). A general formula for the dynamical phase obtained by 
this transformation is given by Eq. (39). In Eq. (40) we have obtained the 
result for im(\l/ in |\l/ ou t) and have shown that for an infinitesimal transfor- 
mation the geometric phase vanishes. We find therefore that the transfor- 
mation given by the matrix Q of Eq. (38) produces a geodesic curve in its 
general meaning. Cyclic transformations for the Pancharatnam phases have 
been obtained by the use of eigenvalues and eigenvectors of the transforma- 
tion matrix but such transformations give only dynamical phases. Geodesic 
'curves' fulfilling the geodesic equation (24) are obtained for special transfor- 
mations. For such geodesic 'curves' the 'horizontal' property of the 'parallel 
transport' is destroyed but not the geodesic property. Possible phase jumps 
of 7r are related to a certain discontinuity in the geometric phase. Although 
the geodesic equation (24) is fulfilled only for special cases, describing the 
development of two levels out of the three-level system, it has been shown 
that the general transformation of the three optical levels by the matrix Q 
is also geodesic in its general sense. Such geodesic 'curves' are obtained for 
vanishing geometrical phases for small changes in the transformation matrix, 
but obtained as a global geometric phase change in the total curve, related 
to non-additivity of the geometrical phase. This geodesic property has been 
justified mathematically also by replacing the geodesic equation (24) by the 
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present geodesic equation (48). In this new equation only the imaginary part 
of the Q matrix is taken into account, since the real part of the Q matrix 
does not contribute to Pancharatnam and geometric phases. The use of the 
general equation (18) for calculating geometrical phases is demonstrated by 
a calculation for a special case. 
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